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INTRODUCTION 


Hl 


Kva  Illation  of  the  expected  performance  of  communications  systems, 
sucli  as  the  World-Wide  Military  Command  and  Control  Systems  (WWMCCS),  in 
the  highly  stressed  cnv’ironments  associated  with  nuclear  warfare  is  an 
important  aspect  of  the  development  and  deployment  of  communications 
systems.  In  addition  to  being  subject  to  the  physical  damage  produced 
liy  nuclear  detonations,  communication  systems  are  additionally  sensitive 
to  degradation  of  the  propagation  medium  through  which  their  signals  are 
transmitted.  .An  extensive  computational  capability  has  been  developed 
in  response  to  the  consequent  need  to  determine  tlie  characteristics  of 
propagation  degradation  and  its  effects  on  communications  systems.  The 
work  reported  here  addresses  one  aspect  of  tliis  degradation  that  had 
seemed,  until  recently,  of  I'elatively  minor  importance.  Recognition  of 
circumstances  under  whicli  propagation  degi'adation  has  gi'cater  significance 
than  formerly  thought  prompted  this  examination  of  its  characteristics. 

Specifically,  we  consider  the  propagation  of  \’ery-low-l  requcncy 
(VI, I,  generally  R to  30  kHz  but  here  limited  to  about  10  to  30  kHz)  and 
1 on  - f reciuency  (I.I',  generally  30  to  300  kHz  but  here  limiti'd  to  about 
30  to  100  kHz)  radio  signals  througli  regions  disturbed  by  t lie  detonation 
of  many  large-yield,  nea  r-sur  1 aci'  nuclear  weapons.  Sucli  eiu-i  ronmiuit  s 
typically  arise  in  the  context  of  large-scale  attack  scenarios  that  also 
include  high-altitude  detonations.  The  sevi'fity  of  \l,l  I.I  propagation 
degradation  associated  with  h igh-a 1 t i t udi'  detonations  generally  is 
su  I I ic  ic  n t 1 y givat  that  careful  coii'- ide  ra  t ion  of  anv  additional  efit'cts 
associated  with  nea  i--sui- 1 ace  detonations  lias  in  the  pa  ^ t been  uniU'Cessary. 
However,  at  SRI  we  hai’e  recently  consideri'd  scenarios  in  v.hich  large 
numlK-rs  ol  1 a rg< ■ -y i e 1 d . near-surface  detonations  occur  in  areas  that  are 
di'Void  of  high-altitude  detonations.  I'urt  hermoi'c' , increasing  VI. I 1.1 
svstem  margins  are  Making  les~  certain  tlv  ciimmun  i ca  t i .uis  ilisruption 
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caused  by  high-altitude  detonations  alone,  which  makes  it  important  to 
include  in  the  analysis  all  significant  cont ributoi’s  to  tlie  propagation 
disturbance . 

The  basic  meclianism  through  which  nuclear  detonations  affect  VLF/LF 
propagation  is  the  same  for  both  high-  and  low-altitude  cietonat ions ; 
ionization  of  the  atmosphere  in  the  vicinity  of  the  detonation.  This 
ionization  is  produced  through  the  interaction  of  the  atmosplicre  with 
energetic  radiations  from  tlie  detonation  and  from  its  dcl)rls.  High- 
altitude  detonations  tend  to  illuminate  the  atmosphere  from  above  Ijecause 
the  atmosphere  about  the  burst  is  too  thin  for  appreciable  interaction 
witli  it  to  occur.  Tlie  i-arliations  from  near-surface  detonations  however, 
tend  to  be  contained  by  the  atmosphere  surrounding  the  detonation.  The 
radiation-caused  ionized  rc'gions  consequently  differ  substantially  in 
form  for  the  two  detonation  types. 

The  ionization  associated  with  a high-altitude  detonation  consists 
of  a relatively  unifoi’m  layei’  «hioh  cxtcnii>^  over  a laj-gt’  ai-ea.  The 
slow  lateral  variation  of  this  layer  permits  treatment  of  its  propagation 
effects  in  terms  of  propagation  beneath  laterally  unifcirm,  oi'  stratilied, 
ionization  layers.  ’fhe  tlu'ory  of  such  propagation  has  becMi  thoi-oughlv 
developed.  It  is  this  approach  that  underlies  presently  a\'ailable  VI. 1 1.1 
propagation  codes. 

i'he  ionized  region  produced  by  a near-surface  detonation  is,  in 
contrast,  initially  confined  closely  about  the  burst.  Iveii  though  the 
ionized  r('gifin  --preads  upward  and  outward  as  the  debris  risi’S,  the  lateral 
variation  of  the  ionization  gem'rally  precludes  an  adecpiati’  treatment  ol 
its  propagation  effects  in  terms  of  propagation  Ixuieath  stratilied  layers, 
■fhe  approach  adopted  in  this  study  to  treat  this  environment  is  t con- 
sider th<’  ambient  ionosphere  and  the  additional  ionization  producefl  by 
th<'  detonations  as  a rough  boundary,  oil  which  the  transmitte<l  signal 
scattei’s.  Ill  ve  1 opiaent  ol  this  modi>l  reciuires  con  s idi  ‘ ra  t i on  of  several 
fac  tors. 


In  Section  1 1 of  this  report  we  examine  the  characteristics  of  the 
ionization  produced  by  near-surface  detonations  more  closely  and  develop 
a statistical  model  of  the  environment  produced  by  many  such  detonations. 
In  this  model,  the  environment  is  characterized  by  the  average  number  of 
detonations  of  a given  type  that  have  occurred  per  unit  area.  In  Section 
III,  we  apply  the  model  to  determine  the  perturbation  to  the  coherently 
reflected  component  of  the  signal  as  a result  of  the  detonations.  In 
Section  IV,  we  examine  the  characteristics  of  the  energy  incoherently 
scattered  by  the  rough  ionospheric  reflector.  Our  conclusions  and  recom- 
mendations are  given  in  Section  V. 


II  propagation  environment 


A.  Physical  Characteristics 

A near-suj-face  nuclear  detonation,  even  of  largo  yield,  initially 
causes  significant  ionization  within  only  a relatively  small  volume. 

The  most  penetrating  radiations  immediately  produced  by  the  detonation 
are  neutrons  and  prompt-gamma  photons,  with  ranges  of  a few  hundred 
of  meters  at  sea  level.  Thus,  the  initially  ionized  region  about  the 
detonation  has  a radius  of,  at  most,  a few  kilometers.  Sucli  regions 
are  too  small  for  a single  detonation  (or  even  many  detonations  separated 
by  distances  appreciably  exceeding  the  size  of  the  ionized  regions)  to 
affect  markedly  \T,F/bE  signals  unless  cither  the  transmitter  or  tlic  re- 
ceiver is  located  directly  adjacent  to,  or  within,  such  a region. 

As  a consequence  of  its  radioactivity,  however,  the  detn’is  from  the 
near-surface  detonation  is  also  a significant  ionization  source.  Tliis 
debris  rises  to  substantial  altitudes  if  the  yield  is  sufficiently  large. 
The  range  of  ionizing  radiation  from  the  debris  generally  is  gi'eater  in 
the  lower-density  air  of  higher  altitudes;  f urt liermore , radiation  directed 
upward  encounters  decreasing  air  density,  and  thus  increasing  range,  as 
it  progresses.  These  factors  suffice  for  the  gamma  radiation  from  (ieiiris 
at  altitudes  aljove  about  20  km  to  produce  an  ionized  column  extending 
upward  from  ttie  det)ris  to  the  amtJient  lower  ionospliere.  Altitudes  of 
this  order  are  attained  by  deliris  from  near-surface  detonations  with 
yields  greater  than  about  1 MT . The  lateral  extent  of  the  dettris  ioniza- 
tion also  grows  with  increasing  altitude  since  the  intensity  of  tlie  ioni- 
zing j'afliation  varies  less  and  less  rapidly  with  latc'ral  distance  as  the 
height  above  tiie  debris  increases. 

Tliis  behavior  is  illustrated  in  figure  1 by  contours  in  a vertical 
cross  section  through  the  detonation  point  of  a constant  electron  density 
of  10  cm  at  various  times  following  a 2-M'f  (k'tonation  at  2.5  km  altitude. 
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OF  10  el  cm-’  DENSITY  AT  VARIOUS  TIMES  FOLLOWING  A LOW-ALTITUDE  DETONATION 


This  ionization  level  provides  a roup;h  indication  ot  a region  sut  i ic  i ent  1 y 
ionize’  to  affect  \T,F/L1'  signals.  Kven  3 minutes  following  the  detonation, 
the  top  of  the  debris  cloud  has  risen  suf  I'icient  ly  (about  1 1-km  altitude) 
to  produce  significant  ionization  above  it.  Tlie  dc'bi’is  continues  to  i-ise, 
with  some  lateral  spreading,  until  about  20  minutes  following  the  rletona- 
tion.  By  this  time  it  has  stabilized,  witli  its  top  at  an  altitude  of 
about  21  km.  The  size  of  the  ionized  regii^  , grows  until  tlie  (kljris  has 
staltilized;  subseciuent  ly , the  size  of  the  ri’gion  containing  significant 
ionization  may  continue  to  grow  as  the  debris  slowly  spreads  laterally. 

The  ionization  intensity  gradually  sul)sides,  however,  a.s  1 ’’.e  detjris 
i-adioact  i vity  diminislies  . 

Calculation  of  the  effects  of  deliris  ionization  on  a signal  invol\-es 
cotis idera t ion  of  both  the  electrons  and  tlie  lieavi-  • positive  and  mgati\-e 
ions  within  the  ionized  region.  A numtHO’  of  simplify’  .g  ;i  ppi-ox  ima  t ions 
a !■('  possitile  for  tlie  altitude  ranges  of  primary  inti'rest  in  this  study. 

,\t  altitudes  in  tlie  atmosphere'  below  about  70  km,  collisions  between  the 
ionized  sp<’Cies  aiul  neutral  atmospheric  molecules  occur  with  sufficient 
rapidity  that  the  propagation  environment  is  well  descrilied  by  an  isotropic 
conductivity,  ;.  I'lie  rc'fractive  index  of  the  medium,  n,  is  given  as  a 
function  of  tlie  conductivity  by 
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( 1 ) 


where  s is  the  permittivitv  of  free  space  and  is  the  wa v<  angular 
o 

1 reciuency  . 


In  forms  of  the  elec-iron  and  ion  number  densities, 
t i ve 1 V , * the  conductivity,  in  turn,  is  given  by 


\ and  \ 

I i 


res  pec  - 


N . = 2N  • \ (S  the  negative-ion  densitv)  is  the  sum  of  the  mg.. five 
,ind  positive  ion  number  densities.  The  negative-  and  positive  ions  are 
assumed  to  have  a pp  rox  i m;i  t e - 1 y ee|ual  masse-s  and  i on-ne-ut  ra  1 collision 
ra  t es  . 
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where  c is  the  maRnitiitle  of  the  electron  charge,  , ami  are,  respec- 

e 1 

tivelv,  the  electron-  and  ion-neutral  collision  I'ates  and  m and  m are 

e 1 

their  respectiv'e  masses.  A useful  conductivity  parameter  is  the  quantity 

= 0/e  . Ci) 

r o 

Tliis  parameter  has  been  used  to  characterize  the  ionization  envii’onment 
in  most  of  the  calculations  to  be  described. 

Reflection  from  the  ionospliere  involves  penetration  of  the  incident 
wave  into  a medium  whose  properties  do  not  clianse  sufficiently  abi'uptly 
to  define  a precise  reflecting  surface.  It  is  possible,  however,  to 
determine  the  region  in  wliicli  reflection  of  the  wave  primarily  occurs. 

This  region  is  characterized,  for  altitudes  below  about  70  km,  by  values 
of  ^ ^ that  appro.ximate  the  wav'e  angular  frequency,  . Koi-  a wave  incident 
at  an  angle  G from  the  vei’tical,  for  example',  field  and  I iigel  find 

c = cos^G  ( 1 ) 

provides  a useful  reference'  surface'  uitliin  this  rt'gion. 

The'  pre'cise  locatie^n  of  such  a sui'face'  is  ne)t  critical  if  tlu'  surface' 
is  properly  used;  we  have'  suppressed  tlu'  angular  de'pendence  in  Kq  . (1) 

to  simplify  matte'i's  wlien  consielering  a range'  of  i nc  ide-nce'  angles.  Ue 
ge'nerally  iiavc  use'el 


(5) 


t ei  define'  a re'fe'j'ence'  sui'face'  at  altitude's  lu'low  atiout  70  km. 
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In  Figure  2,  vertical  cross  sections  tlirouK'h  such  sui-faccs,  as 

■1  -1 

clefinocl  by  Eej . (51  for  a frequency  of  20  kHz  (t  = 1-  > 10  s anti  thus 

c ^ = 2-  / 10 ^s  ^),  are  shown  for  various  times  after  a 2-Ml  detonation 

at  2.5  km-altitude.  The  correspond iiif^  surfaces  of  constant  electron 

density  were  illustrated  in  Fifrure  1.  The  constant -conductivity  surfaces 

can  be  seen  generally  to  have  less  extreme  vai'iation  in  altitude  with 

2 

distance  from  the  detonation  point  than  do  those  for  10  el  'cm  . This 
difference  results  from  the  reduced  contribution  of  the  electrons  to  the 
conductivit,  at  the  lower  altitudes,  where  electron  collision  rate  with 
neutral  molecules  becomes  very  high. 

If  the  curvature  of  the  conductivity  stratifications  within  tlie 
reflection  region  is  sufficiently  small,  the  relation  of  reflected  to 
incident  wave  can  be  described  by  the  reflection  coefficients  for  an 
etpiivalent  plane-stratified  medium.  The  planar  stratifications  of  the 
latter  arc  assumed  to  be  tangent  to  the  curv'ed  stratifications  of  the 
actual  medium  in  the  region  of  interest.  Specifically,  to  use  this 
approximation  it  is  necessary  that 


k(r  ) = 2-(r  ) r-  1 : (6) 

m i n m i n 

that  is,  the  product  of  the  minimum  I'adius  of  curvature  of  the  sui'face, 

( 1-  ) , and  the  wavenumber,  k,  must  exceed  unity. 

‘ ’ m i n 

To  evaluate  Fd| . (b),  we  note  that  tlie  product  of  t lie  two  principal 
(minimum)  radii  of  curvature  of  a sui-fncc-,  for  example  r^  and  r.^,  is 

'c 

given  generally  bv 


•/  Z - /. 

XX  yy  XV 

wlu’rt'  /.  Z(x,v)  is  the  surlacc'  hi'iglit  as  a function  of  the  horizontal 
coordinates  x and  v.  The  suliscripts  denote  partial  d i f I erent  i a t i on  with 
respect  to  the  indi<uitetl  i-ariable  or  \ariables.  In  our  case,  the  heights 
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AT  VARIOUS  TIMES  FOLLOWING  A LOW-ALTITUDE  DETONATION 


of  the  constant -conciuc t ivi ty  strata  depend  only  on  the  liorizontal  distance 
2 2 1^2 

r = (x"  + y ) from  the  detonation  point.  In  terms  of  r,  Kci . (7) 

becomes 


Z 

r r 


(K) 


To  evaluate  r^  and  r,^  separately,  we  note  that  one  ol  tlie  principle 
radii  of  curvature  for  an  ax.imtitha  1 ly  symmetric  surface  will  he  found 
alon^  the  radial  direction.  The  curvature  for  this  direction  is  niven 
simply  by  that  of  a curve  in  one  dimension;  for  Z = Z(r)  we  have 


(9) 


'fhe  other  principle  radius  is  then,  from  i:qs.  (S)  and  (i)). 


(10) 


Table  1 nivt's  value's  of  r^  and  r,^  calculaU'd  from  I.q  s . (9)  and  (10) 

for  t h<'  , 5,  and  10  minute  con  st  a n t -conduct  i v i ty  contours  sliown  in 

Fiuurc'  2.  'fhese  ctjnduc  t i vi  t y surfaces  defiiu’  refere'nce  surfaces  for  a 

wave  freciuency  of  20  kllz  , or  ' = la  km.  Thus,  Kc)  . (0)  is  satisfied  in 

this  case  for  value's  eil  r ^M-e'atei'  tlian  afietut  2.1  km.  A 1 1 lieMif;ti  seime- 

e ' 

scatte'r  re'Sults  in  t tie  value's  ei  1 r^  anei  r,,  in  'fafile'  1 as  a consetiue'iice'  eil 
tlie'  re'lativelv  e'eiarse'  elistane'i'  incre'iients  at  whicli  tlu'  e-einefuc  t i vi  t ie' s 
wefe'  eh  ■ t e' rm  i ne ’ll  , ffie'se'  value's  (^e'lierallv  e'an  tie-  se'e'ii,  w i t li  tfu'  peissilile 
eve-e'pt  iein  eil  tfU'  t fii'ee ' -m  i ’lU  t e ■ ea.se',  to  '-uppeirt  tile'  Use'  eif  planar  I'e- f 1 e-c  t i ein 
e'e>e' f f i e- i ( 'll  t s . Ne>  d i f 1 i e' u 1 t i i 's  are'  ;i  n t i e- i pa  1 e'el  at  hirfur  1 ree|ue'iu' ie' s , teir 
'.etiie-fi  file-  allei-.'i-d  e'urvature'  is  e>\'i'n  itre-ate'r,  lint  some'  e'aution  is  indieati'd 
in  aelopt  i nr  tlie'  use  e)l  planar  ri' I le  e' t i eui  I'oe  ■ I I i e- i e n t ^ to  e'ha  ra  e' t e • r i re  ■ I ou- 
a I t i I uele-de' t rina  t i ein  el  feet'  .it  I r< '9  ue ' ne' i i s fielow  20  kit.'. 
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Table  1 

Principal  Radii  of  Curvature  of  Constant-Conductivity 
Surface  for  20-kHz  Sifjnal 


Distance 

from 

Detonation 

(km) 

Radii  of 

Curvature 

At  3 m i n 

At  5 

min 

At  10  min 

Pi 

(km) 

P2 

(km) 

Pi 

(km ) 

P 2 

(km) 

Pi 

(km) 

P2 

(km) 

2 

2.5 

(5.5 

50. 

07. 

151  . 

35  . 

5 

10.4 

8.3 

14.5 

48. 

000. 

78. 

10 

138. 

40. 

33. 

29. 

--  , 

50. 

20 

580. 

97. 

94  . 

40. 

126. 

57  . 

50 

000. 

320. 

240. 

250. 

250. 

3 1 0 . 

More  complex  expressions  for  the  refractive  index  than  F!ci . (1)  are 
necessary  to  characterize  the  ambient,  or  near-ambient,  nifttittimc  iono- 
sphere, for  which  ionization  levels  are  too  low  below  70-km  altitude  to 
reflect  \'bF/I,I'  signals.  Definition  of  a useful  reference  height  is  also 
less  straightforward.  Tlie  signal's  interaction  with  the  ionospliere  can 
still  be  described  adequately,  liowever,  in  terms  of  planar  reflection 
coefficients,  and  these  coefficients  can  be  readily  calculated  foi-  tliis 
environment  as  well  ns  the  simpler  daytime  and  strongly  disturbed  environ- 
ments to  which  Kqs.  (1)  and  (2)  apply. 

We  have  calculated  electron  and  ion  densities  for  detonations  ix'- 
presentative  of  tliose  of  interest  using  primarily  tiu'  weapons-e  ffect  s 
phenonu'iiology  code  W'KPII  V,  whicli  was  the  most  recent  version  of  tliis 
"ontinually  t'volving  code  av'nilable  to  us  at  the  onset  of  our  study. 

.Some  results  obtained  with  its  successor,  WKPII  VI,  will  also  be  deserilH'd 
briefly.  The.se  latter  results  were'  obtained  in  application  of  the  model 
to  a particular  cMivironment  currently  being  studied  under  another  DN.\ 
contract.*  W'KPH  VI  inclufUs  improved  modt'ls  of  late-time  di'bris  spread 
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font  rac  t DN.AOO  1 -77-C-OOl  H . 


uliich  are  potentially  j'clovant  to  low-a  1 1 i tiidc-rlet  ona  t ion  onvij-onmcnts  . 
Howevci’,  these  improvements  do  not  f>i’oatly  affect  the  debris  beliavior 
diirinff  the  specific  time  period  to  wliich  the  fireatest  attention  has  been 
dii'ected  in  this  study,  namely  the  few  tens  of  minutes  during  which  the 
propagation  effects  arc  expected  to  peak  as  the  debris  stabilizes. 

IJ . Statistical  Models 

Given  the  moans  for  calculating  elocti’on  and  ion  densities  in  tlie 
disturbed  environment,  the  calculation  can  be  very  laborious  if  approaclied 
simply  as  a matter  of  adding  up  tlie  contributions  to  the  total  ionization 
production  that  result  from  each  detonation.  An  alternative  approacli  has 
been  adopted  here  to  avoid  the  task  of  producing  a detailed  description 
of  the  ionization  distribution,  and  consccjuent  propagation  effects, 
caused  by  different  detonation  placements.  Tlie  average  effect  on  propa- 
gation of  an  ensemble  of  randomly  placed  detonations  will  instead  be 
determined;  the  ensemble  is  characterized  by  the  average  density  of 
detonations  per  unit  surface  area.  Althougii  some  clustering  of  detonation 
locations  might  bo  anticipated,  wo  here  address  only  the  simpler  problem 
of  uniformly  random  detonation  placements. 

For  sufficiently  low  detonation  densities,  the  regions  affected  by 
the  indiv'idual  detonations  seldom  overlap.  If  z<'ro  overlap  is  assumed, 
av'erages  over  tlie  ensemble  of  burst  locations  can  be  evaluated  to  good 
approximation  in  terms  of  an  average  over  the  i'<>gion  affected  by  a single 
detonation  and  the  fraction  of  the  total  surfact'  ari'a  that  all  such  region 
const  itut('.  For  the  purpose  of  i 1 1 ust  ra  t i v'e  calculations,  we  have  also 
assumed  that  all  detonations  have  identical  cha rac t er is t ics  and  occur 
simultaneously.  With  these  assumptions,  calculation  of  the  effects  of 
but  a single  isolatt'd  detonation  provi<ies  the  essi'iitial  information 
neci’ssary  to  evaluate  the  desiiad  f|uaMt  i t j<s  . 1 he  limitation  to  a single 

type  of  (U'tonation  can  riadily  be  relaxed,  il  necessary,  by  considering 
a set  of  pi- rt  urb('<l  regions,  one  lor  each  detonation  tvpe. 

fh('  I'cgion  alfected  by  a single  iletonation  can  be  coii'-idered  to  have 
a linitc’  a I'ea  , even  though  the  rig  ion  is  not  delimited  by  a '-udden,  -harp 


cutoff  of  ionization  production.  The  final  result  is  not  sensitive  to 
tile  specific  boundary  assumed  for  this  region  as  long  as  the  region  is 
large  enough  to  encompass  the  significant  cietonation  effects.  The  prob- 
ability density,  p(x),  tor  some  quantity  of  interest,  x,  can  then  be 
expressed  as  a weighted  sum  of  the  density  for  the  affected  region, 

p,  (x)  and  that  for  the  unaffected  area,  p (x); 

1)  o 

p(x)dx  = [p.\  ‘ (1  - p.A)p^(x)_dx  , (11) 

where  p is  the  density  of  detonations  per  unit  surface  area,  and  .)  is  Uie 
area  of  the  region  taken  to  be  affected  by  the  detonation.  It  is  evident 
from  Kq . (11)  that  the  product  p.A  must  be  much  less  than  unity  il  this 
approach  is  to  be  meaningful,  and  the  contribution  of  to  p(x)  will 

generally  be  small  in  such  cases. 

Averages  are  similarly  constituted: 


/> 


xp  ( X ) dx 


(x)dx 


(1 


(x)flx 


(12) 


The  probability  density  for  the  ambient  envi  I'onment  , p^^(x),  will  g<‘nei-all\- 

be  a delta  function  since  the  quantity  x will  have  some  specific  \'alue, 

X , for  tlie  undisturbed  environment.  I'hus,  t lie  rigiit  integral  in  K(|  . (12) 

o 

can  be  evaluated  to  give 


I 


X = xp  (x)dx  * (1  - pA)x 

■ |)  o 


(12) 


■fliis  general  aiiproach  can  be  ext('nded  to  <let<5nation  densities  gi-eatei- 
tlian  those  to  wliich  i:<|s.  (11)  and  (12.)  can  1m'  applied,  and  for  wliich  t tie 

disturbed  efiects  can  indeed  liecomi'  significant.  I.\'cn  wlien  the  regions 
alfected  Ijy  the  dillerint  detonations  overlap  a ppn  ■<■  i all  1 y , t lie  area  in 
which  two  or  more  chdonations  contribute  similarly  to  ttu'  total  ionization 
often  ri'inains  only  a fixiction  of  the  total  area  ol  ov<‘ila|i.  In  the  re- 
miiimh'f  ol  the  overIa|niing  region,  the  ionizing  flux  from  one  ol  the 
detonations  is  markc’div  more  intense  than  that  from  the  otlno's,  and 


detonations  is  markc’div  more 


consideration  of  only  this  dominant  ionization  source  should  provide-  a 
rea sonable  approxima t ion . 


Although  we  have  not  evaluated  the  accuracy  of  this  approximation 
quantitatively,  it  appears  to  be  useful  to  liifjher  eletonation  dc-nsities 
than  mifiht  initially  be  thout^ht  to  t)c  the  case.  Consider  I h<-  liorizontal 
e-foldinR  distance  for  the  ionizing  flux  from  a detonation,  v.liich  is 
shorter  at  lower  altitudes  tlian  at  hip;hei-  altitudes.  i he  altitude  raiiRc- 
of  interest,  and  correspondingly  tlie  e-foldins^  distance  of  consequence, 
decreases  with  decreasing  distance  from  a detonation.  The  averatie  linear 
distance  ov'er  which  dominance  in  ion  production  shifts  from  one  de-tonation 
to  another,  in  the  region  midway  between  two  detonations,  consequently 
decreases  with  increasing  detonation  density.  It  is  tills  intermediate 
region  that  consideration  of  only  the  dominant  detonation  liadly  mis- 
treats. Thus,  the  hardening"  of  the  boundaries  between  dc-tonations 
with  decreasing  separation  betwc'en  tliem,  in  reducing  the  area  associated 
with  the  boundary  region,  serves  to  extend  the  applicability  ol  t he 
approximation  to  higher  fletonation  densities  tlian  would  otherwise  have 
been  the  case. 

An  efficient  appi’oach  to  implementation  of  thi-  motk'l  is  to  deter- 
mine the  probability  distribution  for  the  horizontal  distance  to  the 
nearest  detonation.  Definition  of  the  functional  dependencf'  upon  this 
rlistance  for  any  quantity  of  interest  then  ]<'arls  immediately  to  its 
probabilitv  distribution.  ihe  distance  distribution  function  derives 
from  the  Poisson  distribution,  which  describc'S  the  probability  for  t he 
number  of  independent  events  within  a region,  gi\cn  their  a\-erag(-  density. 
To  obtain  the  distribution,  we  note  tluit,  for  small  t'nough  r<gions,  t he 
J probability  that  theri-  will  be  one  event  within  the  rt’gion  bi-comes  t he 

f densitv  ol  events  (number  per  unit  a n-a  ) , P,  t inu'S  the  area  ol  t he 

, region.  A; 

I 

j Ptevent  in  A)  P(A)  pA  . (Ml 

' Converselv,  the  probabilitv  that  t he  I'e  will  be  no  i vents  i ji  A (s 

ihe  probabilitv  that  there  be  more  than  one  event  in  \ bei'omes  r.igligibb 
' relative  to  the  probabilitv  o(  one  event  in  A (or  A s u ( 1 i e i eti  t 1 v t::a  1 1 . 


Ih 


P(no  event  in  A)  = P(A) 


1 


(15) 


- P-\ 


For  areas  of  appreciable  size,  we  can  write 


N 

P(A)  = lim  " P(A  N)  (Ki) 

N 


that  is,  we  divide  A into  N subroRions  and  consider  the  limit  as  N becomes 
lar^e.  In  this  limit,  however,  the  snbref^ions  l)ccome  small,  and  wc  can 
employ  Eq . (15)  to  determine  P(A/N).  Thus,  takinp;  P(A/N)  - 1 - pA/N,  we 
have,  from  Eq . (16), 


N 

P(A)  = lim  ",  (1  - . (17) 

Equation  (17)  is  merely  one  way  of  defining  tlie  exponential  function, 
and  wc  find 


P(A)  = e , (IS) 

the  Poisson  dist  rilnit  ion  for  no  events  in  A.  For  a cii'cular  i-egion  of 
radius  r,  the  prot)al) i 1 i t y of  no  detonations  witliin  the  I’cgion  is  ttuis 

2 

P ( r ) --  e”  . (19) 

Now,  we  want  to  calculate  the  probability  that  tlu'  lU'arest  burst  is 
at  some  distance  r.  For  this  situation,  we  I’equire  both  that  no  (h'tonations 
occur  at  shorter  tiistances  and  that  one'  occurs  at  i-.  Thc'se'  are  independent 
e'V'ents,  so  the  joint  probability  is  t he>  product  e><"  the'  re’spective’  preib- 
abilitie'S.  Thus,  we<  have',  from  Fe|s.  (II)  and  (19),  the  probability  density 
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f line t ion 


p(i-)(ir  = 


‘j-  ^ 


r e 


dr 


(20) 


that  the  nearest  detonation  to  a point  is  a distance  r from  it. 

Kqiiation  (20)  reduces  to  the  expression  for  isolated  detonations, 

Eq . (11),  if  the  density  p is  sufficiently  small.  We  note  that  p(r)dr 
can  fienerally  be  written 

p(r)dr  = ' p(r^.A)clr  , (21) 

where  p(r^.\)dr  is  the  density  function  for  the  occurrence  of  r in  A (and 
p(i'jA)  the  density  function  for  r not  in  A),  and  A is  now  tlie  area  affected 
l)y  a detonation.  Equation  (21)  can  equivalently  ()e  written,  introducing: 
the  conditional  protiabi  1 it ies  p(r|A)  and  p(r|A)  that  r is  within  or  with- 
out A respectively,  as 


p(r)dr  = [p(rl\)P(A)  • p(r|A)P(A)  dr 


But 


P(A) 


whicdi  becomes,  lor  1, 

P(A)  ==  ; = pA 


(22) 


(23) 


(2  1) 


Ihus,  notinfp  P(A)  = 1 - P(A), 


p(r)drs  [jAplrlA)  (1  - ,A)  p(r|A)  dr 


( 2 a 1 


IM 


Ill  COHEliKNT  SIGNAL 


The  roles  of  the  (ideally)  smooth  earth's  surface  and  lower  iono- 
sphere in  the  propagation  of  Ll'/VLK  sif^nals  around  the  earth  tu-tween  t liem 
can  be  conveniently  expressed  in  terms  of  tlieir  plane-wave  reflc'ction 
coefficients.  Tliese  coefficients  Kive  the  ratio  of  the  (ielti  rellected 
1)V  the  boundary  to  that  incident  upon  it.  Decomposition  of  the  si'^nal 
radiated  by  the  source  into  an  angular  spectrum  of  pl.''.nc  waves,  followed 
by  modification  of  cacli  plane-wave  component  by  the  appropriate  j-eflection 
coe  f f ic  ient  ( s ) , yields  the  field  integrals  that  are  evaluated  by  vai'ious 

C 

presently  available  computer  codes. 

The  specular,  or  coherent,  reflection  from  a rougli  lioundary,  such 
as  the  perturbed  ionospheres  considered  here,  can  also  bi'  dc'seribed  in 
terms  of  such  i-ofU'Ction  coefficients,  which,  however,  now  depend  on  t lie 
average  properties  of  the  surface.  An  incoherent  return  is  also  built 
up  as  a result  of  scattering  of  t'nergy  out  <if  the-  specular  direction  by 
the  surface  i i-regul  a ri  t i os . Ibis  scattering,  which  is  a function  of  t he 
variability  of  the  sui'face  pi'operties  about  tlu'ir  average  values,  will 
bc‘  discussetl  further  in  the  next  section. 

It  was  .shown  in  .Se<’1  ion  I 1 -A  above,  that  the  perturbed  retleeting 
surfaces  produced  by  debris  ionixation  I rom  the  large-yiild,  low-altitude 
detonations  of  inteiwst  here  have  su  I f i c- i en  t 1 y large  radii  ot  cuiwature 
that  the  interaction  with  t hf'lii  of  th('  inciihnit  field  can  also  b(‘  described 
locally  in  terms  of  planar  rrllection  coe  I I ic  ient  s . I'his  result  permits 
substtintial  simplification  of  the  calculation  ot  the  ;i\erage  reflection 

coelficients  for  the  pertvirbed  surlac< as  ;i  fir--t  approximation  we  need 

onlv  average  thi'  local  reflection  coellicient  ovir  the  entire  surlace. 

A more'  accurate  evaluation,  whii'h,  how<  vi-r,  will  not  be  tiiken  up  lu  ta  , 
would  include'  me-eeli  I ic;i  t ient  eel  the'  i ne- ieie  n t - t ie' 1 el  stre'tigth  ( I reem  that  in 
the'  unpe' r t u rbeel  e 'iiv  i retiime 'ti  t 1 tee  acceeunt  leer  the'  pee--sibly  signilictint 
I'llie'l-  ell  mtlltiple'  -e'titter  and  shaeleewing. 
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Let  h'  be  t lie  hoiKht  at  which  c satisfies  Kq . (5),  and  li  the  refer- 

r 

ence  heif^ht  chosen  for  the  averafrc  reflection  coefficient,  K(ti,h),  of 
the  roiigh  surface  for  a plane  wave  incident  at  an^le  6 on  the  mean  surface 
Then  the  averaf;e  reflection  coef  jient  is  '5V 


/' 


R(9,h)  = jp(r)H  6,h'(r);r^  exp 


I -1  - cos9  - 
i : h - 


h' (r)  ! dr 


(2b) 


whore  R^[G,h'(r);r  is  the  local  reflection  coefficient,  referenced  to 
height  h',  and  p(r)  is  p;iven  by  Kq . (20).  The  exponential  function  shilts 

the  phase  of  from  its  local  reference  lieiRlit  h'  to  the  heif;ht  li , which 
can  appropriately  be  taken  to  be  the  average  of  tlie  h^  over  r. 

The  ionosphere  can  bo  separated  for  low  detonation  densities  into 
an  unperturbed  ambient  background  and  nonoverlapping  perturbed  rc'gions 
about  each  detonation,  as  described  in  Section  11-13,  above.  'I'hc  dependenci 
of  h'  on  r within  the  perturbed  region  is  determined  by  the  sliape  of  the 
constant -conductivity  surface  for  the  value  of  c given  by  l.ci . (5).  In 

the  ambient  region,  h^  is  a constant.  Thus,  Kq  . (2(3)  becomes 


R ( e , h ) 


exp(  - 


•p  r 


')  R [ c , h ' ( r ) ; r ; exp 
o ■ 


■1-cosG.  , 
1 -7 h 


h'(r)  jdr 


exp(  ■ 


R^(",li.^>  exp 


, 1 cos- 
1 ^ ( h 


(27) 


where  A is  the  area  of  t lie  perturbed  region  and  R^(G,li^)  is  the  ambifuit 

value  of  R , referf'ticed  to  the  ambiei  t height,  h,. 

o A 

\ number  of  calculations  have  been  made  that  demonstrate  the  signi- 

lieanci-  of  the  detonation  effects.  fable  2 presents  the  ri'sult  of 

ivaluating  Kci  . (27)  for  t li('  20-kll/  rc'flection  coelficieiits  ol  a \'ertically 

polarixed  waves  at  vai’ious  times  following  numbers  of  2.a-MT  detonations 

at  2-km  altitude  in  an  ambifuit  da\ t ime  environment  . Results  are  given 

— (1  — a — 1 

for  average  detonation  diuisities  of  i)  (ambient), 10  , It'  , and  in  per 

■-(luare  kilometer.  ’ITie  nature  of  the  i otios|iheric  perturbation  caused  by 
each  detonation  is  illustratefi  by  the  Cfuis  t a n t -e  1 ec  1 ron-detis  i t y surfaces 
plotted  in  figure  1 and  the  c on,- t an  t -conduc  t i i t y surlaces  plottid  in  figut 


Tabic  2 


■Av'cragc  Hof  lection  Coefficients  for  20-kHz  Vertically 
Polarized  Plane  Wave  Incident  at  Anfilc  6.  in  Disturbed 
Daytime  Dnvii’onment 


Detonat ion 

He  fleet  ion 

Cocf  f icient 

COS  - 

1 

Density 

ot  5 minutes 

at  10 

minutes 

at  20  minutes  | 

(kin”“^ 

iKl 

ai'K  K 
(dc'K) 

|h1 

a I'K  H 

(defi) 

|H| 

ars  K 
(dc-K) 

0.1 

0 

0.(;75 

-KiO. 

0 . 6 76 

-160. 

0 . 676 

-160. 

-(S 

10 

0.671 

-169. 

0.678 

-169. 

0 . 676 

-168. 

io“'’ 

0.669 

-1  19. 

0.687 

-1  19. 

0.697 

-1  18. 

10-' 

0.671 

-117. 

0.686 

-118. 

0.671 

-120. 

0.2 

0 

0.  163 

-M  l . 

0.  153 

-Ml. 

0 . 163 

-M-1  . 

-0 

10 

0.1  19 

-Ml  . 

0.  163 

-Ml  . 

0.  163 

-110. 

1()-'’ 

0.1 31 

-123. 

0 . 168 

-121  . 

0.171 

-119 

lo"  ' 

0.333 

-05  . 

0 . 3.68 

-66 . 

0.135 



The  local  reflection  coefficients  were  determined  foi'  this  calculation 
bv  fitting:  CNponential  firadii'uts  at  the  reft'rc'nce  heif;ht  to  the  conductivity 
prol  iles  calculated  for  difi'erent  distances  from  tli(>  detonation  ))oint  . 
Keflection  coellicients  for  exponential  conrluct i v j t y profiles  are  available 
that  (Micompass  the  rant;e  of  conduct  iv'i  I y |^i-adients  found  for  thest'  prn- 
files.  These  I’esults  were  list'd  in  a look-up  table  witli  i n ti' rpo  1 n t i on  on 
pliase  and  amplitude  tietwet-n  tabulaied  values  to  i^t'uerate  the  intcftrand 
of  1 (| . (27).  Tlie  two  incidence  anttles  lOr  wliich  calculations  were  made, 

H 1 . ,’f  (cos  = d.l)  and  7H.r>  (cos  = 0.2),  are  representative  ol  t lie  ramte 
ol  aufilt's  !-iKnilicant  in  moderal  e-t  o-lonp;  distance  VI.F/l.K  propauation. 

fht'  ambient  reflection  height,  rather  than  an  avt'raue  hi'ii;ht,  was 
used  as  the  phase  refert'tice  for  the  d i s t urbed -en  vi  ronmeii  t reflection 
coe  I I icienl  uiven  in  lable  2.  With  this  choice,  any  v-han^te  in  reflect  ion- 
coeflicit'iit  phase  I rom  its  ambit-nt  value  is  not  masked  by  s imu  1 ( ane<nis 
chaiutes  in  1 hi'  relereiice  heiirhl.  \n  average  reference  heiitht,  liowcver. 


would  be  pro  f'ci'iib  lo  in  nny  fui'tber  c;i  Icul  :i  t ions  involving  these  reilec- 
tion  cocl ficients . 

Aside  irom  sf>nie  jitter  in  the  results  whicli  is  caused  primarily  by 

the  coarst'uess  oi  the  integration  stop,  the  data  pi’osented  in  Table  2 

-G  —2 

show  a consistent  patteini.  Detonation  densitic-s  below  10  km  produce 

-o  -2 

litt’e,  if  tiny,  effect.  At  a density  of  10  km  , the  detonations  caus<' 
a perceptible  phase  sliift  and  iwwivict  ion  in  reflection  amplitude  at  f>  min. 
The  phase  shift  reimiins  essentially  constant  Itetween  o ;ind  20  min,  but 
the  timplitude  per  t u rltti  t i on  chant;es  from  the  I’educt  ion  at  5 min  to  an 
enluinc(‘ment  at  20  min. 


Ihe  re  f I ec  t ion-coe  f f i c ien  t perturbation  shown  by  I'able  2 for  a 
- I -2 

detonation  den>it\  ol  10  km  is  substantial.  ,\s  with  that  tor  a 
dfuisity  o!  10  ’ km  , the  ptiase  shilt  remains  nearly  constant  between 
.')  and  20  min.  1 lie  I'el  lection  amplitude  is  markedly  i-educed  relatiye  to 
amltient  at  a min,  with  the  rt'duction  decreasinjt  with  increasing  time, 
lor  cos  , 0.1,  the  amplitude  appears  to  haw  recoyered  neai’ly  t . its 

ambient  yalue  ol  20  min.  llowe\-ei',  the  re  1 lee  t ion-coe  f f ic  ien  t behayior 
for  a detonation  density  ol  lo  km  , coupled  with  the  persistence  of 
t lu'  phase  perturbation,  makes  such  interpretation  suspect;  rather,  the 
amplitude  perturbation  may  well  eyolve  into  an  enhancement  at  later 
times  befori  ultimately  reco\'e  r i nu'  to  its  ambient  yalue.  At  moderate 
(f)i'der  of  1000  km)  distances,  these  chany;('s  in  re  I U-c  t i on-coe  1 1 ic  i en  t 
am|)iituile  translate  directly  into  changes  in  I'ecciyed  lield  strength. 


I xamination  of  the  rel  lection  coefiicients  civen  in  Table  2 as  ,i 
lunction  ol  cos  , the  cosine  of  the  incidence  anele,  shows  them  to  be 
reasonably  well  lit  by  the  function  form 


i;  - ex  p ( ..  c-os  ) , ( :'S  ) 

where  ■ i ^ is  complex,  as  illustrated  in  liyvire  A lor  the  20-min 

data.  ihis  form  has  been  fouiul  to  describi  well  the  re  I 1 ec t i on-coe  I I i c i i 
rallies  at  trra/int;  incidence  anyles  (cos  1)  loi'  a wide  \ariet\’  ol 

i onos  pile  r<  s . Its  suitability  here  makes  immediately  releyant  the  Vl.l 


tCTION  COEFFICIENT  ON  INCIDENCE  ANGLE 


modal  propagation  clia  ractori  st  ics  calculated  t)y  Wait  and  Si^ies,  which 
includes  parameterization  of  the  results  in  terms  of  their  dependence 
on  0 . 

Taljle  ;i  giv'es  tlie  values  of  0 derived  foi’  the  ri'flection  coefficients 
presentf'd  in  Table  2.  Tliese  Q,  wfiicli  arc  an  average  tjf  t lie  values  deter- 
mined separately  from  the  reflection  coefficients  for  cos  v = 0.1  and 
cos  9 = 0.2  in  Table  2,  yield  the  linear  curves  shown  in  figure  2.  The 
values  of  are  referenced  to  a height  of  70  km;  this  height,  although 
somewhat  high,  is  the  lowest  phase  reference  height  for  which  modal  con- 
stants  can  readily  bc'  inferred  from  Wait  and  Spies  data.  Pei-usal  ol 

Wait  and  Spies  results  show  that  tlie  values  of  s obtained  for  detonation 
-a  -2 

ciensities  up  to  10  km  all  correspond  to  attenuation  rates  over  sea 

water  within  0.1  dB  of  2.1  dB/Mm,  the  ambient  value.  lor  a detonation 
-1  -2 

(h'lisity  of  10  km  , the  attenuation  rati'  at  .T  min  is  estimated  to  be 

greati'i’  than  twice  tlie  ambient  v'alue,  with  some  rt'covery,  to  about  twice 

the  ambient  v’alue,  at  20  min  following  the  di't  ona  t ions  . The  pha  si’- ve  loc  i t \ 

departures  from  their  ambient  value  are  significant  foi'  detonation  densitii 
-o  -2 

of  10  km  or  greater. 


Taljle  2 

I’a  rame  t er  i za  t i on  of  K lor  20-kHz  Wave  in  tlie  lorm 

V 

K - - expfci  cos  -.) 

V'  1 


Det  ona  t i on 

Density 

at  5 minuti'S 

at  10  minutes 
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Cl 

Cl 

a- 

Cl 
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I 

H 

1 

i< 

1 

0 

-2 . 9 1 

9.2.5 

-.2.9  1 

<).2;5 

-2.9  1 

9 . 95 

10 

-2.97 

9.  ,5  7 

-2.92 

9.. 57 

-2.<1  1 

9 . 70 

-5 

10 

-1  . 10 

11.22 

-.2  . ,H2 

1 1 . 1 

-:i . 09 
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- 1 

10 

- 5 . f)  5 

10.  .5.5 

-.5.  1 ,S 

10.12 

- I .o.s 

10.21 

+ -1-2 

Thi'  o values  foi-  10  km  fall  well  out'-idi  the  range  ol  tlie  l igui-es 
prepared  by  Wait  and  Spii’S,  and  the  necessary  ext  r.i  po  1 a t i ons  i-ould  be 
performed  only  roughlv. 
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lU'sulls  liavc  also  l)ccn  obtained  for  propa^jation  in  a ni^bt  environ- 
ment, as  part  of  a study  recently  completed  undei-  anotliei'  DNA  contract.* 
Since  no  ni^bt  environments  were  considerc'd  in  tlie  work  performed  directly 
under  tbe  contract  i-eported  iu'i'c,  some  results  fi-om  this  otber  study  will 
be  (lescril)ed  briefly. 

NiU'bl  aml)ient  ionosplieres  do  not  permit  many  of  tbe  simplifications, 
described  in  Section  [ 1 -A , tbat  are  possil)le  for  ttie  more  intense  daytime 
ionospb(>re' . Tliis  greater  complexity  makes  diiwct  numerical  calculation 
of  tlie  ionosplieric  reflection  coefficients,  usinp;  a suitable  computer 
coele’ , pertiaps  the  most  st  r;i  in'lit  forwa  rel  appretacli  tei  their  ele  t ermi  na  t ion  . 

A vt'rsietn  oi'  tbe  procedure'  de've'l  etpt'd  foi'  such  ca  Iculat  ietns  by  Pitte'way 
was  adopted  in  tliis  work.  'file  compute)'  code'  imple'inent  i ntf  I'itteway's 
preice'dure'  was  fed  profile'  data  obtainc'el  with  tVKPll  VI,  tlui  most  re'Ct'iU 
ol  the'  U'Kl’ll  series  ol'  nuclear  pht'nomt'nole)p;y  code's.  'file'  detonation  ebar- 
acte'r  i st  ics  assumeel  in  this  calculation  we're  the'  same  as  those  adopti'd 
lor  tilt'  day  t ime -env  i ronmen  t calculation.  'flit'  variation  of  re'flcction 
characteristics  with  bori/eintal  distance  from  such  an  isolated  dt'tonation 
in  the  nin'ht  environment  is  ext'iiipl  i f ied  by  the'  etata  niven  in  fable  1. 

'file  entries  ttive  the  mattnilude,  at  10  min  followinit  such  a detonation, 
of  the  vertically  polari'/ed  rel  lection  coe  f f ic  it'iit  , U , tor  a l.o-kll/ 
wave  incide'iit  at  an  ann'le  of  SO  . A phase  refert'nce  heietht  c'host'ii  to 
lie  in  the  reflection  rettion  is  also  indicated  in  the'  table.  I'he  relCr- 
eiiee-ht' i ech  t criterion  used  to  determiiU'  this  height  reduces  to  1 ()  . (1), 

at  low  altitudes.  Inclusion  of  the  i nc  i<k>  nce-a  n^' 1 e lactor  of  l.ci  . ( I) 

in  the  determination  of  the  reference  lu'inht  leads  to  much  lower  lu'itchts 
beiU'.ith  the  detonation  debris  than  weri'  obtaiiU'd  for  t lit'  davtime  cast', 
di'-fussed  above,  usitu;  t lit'  simplified  formula  of  iti  . (lit  . 

lOr  the  niftlit  ambient  i tinosphe  re  , the  reflection  coefficient  of 
the  hori/ontally  polari/t'd  wavt',  , and  the  conversion  ctie  I 1 ic  ieiit  s 

i;^  and  It  must  also  be  known  to  iletermine  the  propagation  character- 
istic^; these  '-lUantities  were  all  calculated  and  averaged  over  the  dis- 
turbed environment  to  obtain  the  necessary  data.  This  calculatitin  was 
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Tabic  5 


Attenuation  Rates  tor  ir>-kli/  Uavc  Rropa^a  t i nr: 
Over  UiiKi  in  Disturbed  Nisbi  Knvironment  at 
10  Minutes  I'ollowint;  Detonations 


Detonation 

A t tenua  t ion 

Density 

Kate 

(km  ~ ) 

(dli/Mm) 

-5 

10 

D.  70 

■>  / lo"’’ 

■1  . HO 

o / 10“‘’ 

(t . HI 

IV 


INCOHERf-NT  SIGNAL 


A major  effect  of  surface  roughness  on  a sinual  reflected  from  the 
surface  is  tlie  reduction  in  amplitude  of  the  coherent,  specular  return 
that  was  discussed  in  the  preceding  section.  As  is  evident  from  ttie 
formulation  presented  there,  this  amplitude  reduction  results  in  large 
part  from  phase  interference  between  the  contributions  to  the  I'cflected 
wave  from  various  parts  of  the  surface.  The  energy  lost  from  tlie  specular 
return  in  this  manner  is  scattered  into  ottier  directions  rather  than  b(>ing 
absorbed  by  tlie  boundary.  Conseciuent  ly , it  is  of  interest  to  determine 
the  incoherent  power  scattc'ied  to  the  receiver  by  the  surface  as  well  as 
the  Btri'ngth  of  the  coherent  return.  For  moderate  patli  lengths  (up  to 
about  2000  km),  multiple  rel lections  between  earth  and  ionosphere  play 
a secondary  role  in  signal  propagation.  Fnder  these  conditions  the  re- 
ceived signal  can  be  found  relatively  st  ra  iglit  forwa  rd  1 y by  direct  inte- 
gration of  the  contribution  from  each  element  of  the  relli'Cting  surface. 

We  assume  lurther  that  (1)  multiple  scatter  betweiui  points  on  the 
surface'  can  be  neglected,  ( ’)  the  surlace  curvature'  is  not  e'xce'ssive' 

(as  spe'Cified  by  1 C|  . (0)),  anel  (3)  path-U'iigth  el  i 1' ft'rences  fear  wave's 

scatte're'el  f reim  diffe’reuit  part-  of  the’  surface'  are  sufficiently  great 
to  randeimi/e'  tlu'ir  phase's.  Inele'r  these’  ceineli  t ieins  , inte’gration  over 
the’  re’fle’cting  surface’s  e-an  be’  re’eluce'd  to  a summatiein  eif  the'  peiwe’r 
scatte'reel  by  spe’cular  peiints  on  the’  surface’.  A statistical  average’ 
ove’i-all  the’  possible’  re’lle’Cting  surfaces  tlu’ii  give  s the'  incohere'nt 
scatte'ri’el  peiwer  as  the'  proeluct  o<  • lu’  ave'rage-  numbe’r  ei  f suitably- 
eii-ie-nte’'!  spe’cular  points  peT  unit  surface'  area,  the'  aveiuige'  scattering 
oreiss  se’Ctiein  eil  e'ae'h  specular  peiint,  anel  the’  incide-nt  penee'r  per  unit 
surltice’  are’U.  'fhe'  ave-rage'  scatte’ring  creiss  se’e'tiem  can  the-n  be’  e'X- 
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'I  lie  specular  points  are  determined  by  the  slope  of  the  roupli  re- 
tlectinst  surface;  i.c.,  tlic  disturbed  ionosphere.  The  statistics  of 
til  is  surface  are  used  to  calculate  the  av'crane  number  of  specular  points 
per  unit  area  and  the  averasc'd  product  of  the  principle  i-adii  at  the 
specular  point.  In  our  statistical  model  (Section  II-B),  those  cliaracter- 
istics  are  determined  by  those  of  an  isolated  perturbation,  which  can 
b('  evaluated  to  whatever  precision  is  required.  Ue  have  examined  in 
this  initial  effort  only  the  simplest  case  of  well-separated  detonations, 
for  which  the  isolated-detonation  characteristics  are  related  to  those 
of  the  entire  surface  by  Kc; . (11).  1 he  incoherent  received  power  in  this 

case  is  expected,  and  indeed  is  found,  to  be  small  relativ'e  to  tliat  in 
the  specular  return. 

Probability  distributions  for  the  reflectiiiK  surface  slopes  arc 
calculated  directly  from  elt'ctron  and  ion  contours.  Surface  character- 
istics liave  been  dett'rmined  for  3,  5,  10,  and  20  minutes  after  the  detona- 
tions occvir.  Two  options  for  fk'finintt  the  reflectin;^  surface  are  available 
in  the  calculation.  ’fhe  first,  and  poorer,  approximation  is  to  use  the 
lieiq^ht  of  (‘lectron  dtuisity  ;it  10  el  cm  as  a function  of  location.  i'he 
second  option,  wliich  pro\'id<‘S  a mori'  accurate  description  of  the  surface, 
is  to  use  a constant  conductivity,  the  value  of  whicli  depends  on  tlie 
input  lfe(|uency,  as  ttiven  by  Iti,  (5). 

S p<’c  i f ica  1 1 y , the  rms  strength  of  tlu'  vertical  entmponent  , I , <■)! 
the  incoherent  field  is  ejven  by 


v •) 


V.  he  re 

1,^  ineidi'iit  lield  strenr:th  at  the  scattering;  surlace. 

j --catteriny  cross  sectioTi  per  unit  surface  arini. 

r’^  distance  I rom  the  scattering'  element  to  the  reception 

po  j n t . 


in 


nn^lc'  between  tlic  v'erticnl  nnd  Die  incident  ray  at 


the  receiver.  Because  planar  nx‘C)metry  is  used,  this 
aiifile  is  also  the  projection  on  a vertical  plane 
thi’OUKh  the  receiver  and  the  reflection  point  of  the 
scattered-ray  alible  at  the  reflecting  surface. 
d.A  = scattering  surface  increment.  'I'lie  integration  is 
perfonned  over  tlie  scatteritig  surface  .A. 

The  incident  field  at  the  scattering  surface'  is 

V s i n 9 


wlif're' 

1-  = distance  from  the  transmitte-r  to  tlic  scattering  ('lenient. 

V = source  strength,  expressed  as  a voltage. 

0 

sin  V = pattc'rn  factor  fi^r  radiated  field  assuming  a vertical- 
i 

dipole  source,  where 

s = angle  tx'tween  tlu'  vertical  and  the  rav  bf'ing  transmitted 

1 

Because'  planar  geometry  is  used,  tliis  angle  is  also  t lic' 
projection  on  a vertical  plane  threxigh  t lie  transmitter 
and  tile  re'l  lection  point  ol  the  i nc  ideiit -ray  angle  at 
the  reflecting  surface. 


■file  surface  integration  is  accom))  1 i shed  in  polar  coordinates.  lo 
illustrative  purposes,  a planar  propagation  geometry  has  lieeii  ehoseii, 
in  which  Ki.  ( id) ) becomes  (also  using  Ke|  . (dll)) 


I J 


-in"-  . ( ) 

1 ■ 


/ 


( 


wh( ' ri  ■ 

,,  - distance  I rom  the  t ra  iisin  i t 1 1 r to  the  ijoint  ol  relliction 

measured  on  the  ground  plain  . 
h average  height  ei|  the  I'e  I lect  i ng  siirlaee  . 


Barrlck  derives  tlic  average  scattered  power  as  a iunction  of  the  average 
number  of  specular  points  on  the  reflecting  surface  multiplied  by  tlie 
average  curvature  at  these  points.  He  shows  that  this  quantity  is  equal 
to  the  probability  density  for  tlie  slopes  necessary  to  produce  a specular 
point  multiplied  by  an  angle  factor; 


whe  re 


» <- 

X ''V 


angle  between  the  mean  normal  to  tlie  surface  and  the 
local  normal  to  the  surface, 
height  of  the  reflecting  surface. 

surface  slope'  at  the  specular  point;  de'r  i vat  ive's  are 

in  the  .x  and  v directions  respect  ivelv . (1  and  I 

“X  ' y 

arc  computed  in  a coordinate^  system  where  deviations 
arc  measured  from  tlie  average  reflc'Cting  surface  he'ight.) 
probability  of  occurrence  on  tlie  rcfle’cting  surface  of 
a point  having  the  slope  • 


Assumptions  about  the  perturbed  surface’  that  allow  for  the’  analytic 
simplification  are  the’  following: 

(1)  The’  burst  pe'rturbat  ions  arc  symme’tric. 

(2)  All  bursts  are  iek'iitical  and  occur  at  once-  so  that  all  pe-r- 
turlie-d  are’as  are-  ide-ntical  at  any  spe-cified  time. 

(,'l)  The’  e-ffe-ct  of  a single’  burst  e-xte-nds  ove’r  a finite-  are-a. 

■fhe  pe-rturbe-el  re’gion  has  be-e-n  ele-fine-ei  to  e’xte-nel  out  to 
the-  peiint  whe- re-  the-  slope-  be-conie’s  le-ss  than  0.02.  (The- 
ambie-nt  re- 1 1 e-e- 1 i ng  surface-  has  a /e-ro  slope-  in  this  model.) 

(1)  The’  burst  leicatieiiis  are-  re-pre-se’nt  e-d  by  a unifoiTi  raneieim 
spatial  el  i st  r ibut  ion  . 

(.A)  No  ove-rlapping  bursts  are’  alleiwcel.  (The-  probability  of 
this  occur i-e-nce-  is  compute’el  as  a che-ck.) 


If  the-  shape-  eif  the-  pe’ I't  urbe-el  re-flecting  surface-  is  kneiwn,  (i.e-., 

' anel  1!  as  a Iunction  of  x anel  v)  the-  preibabilitv  ei  i st  ribut  iein 
• X ' y 

p('  ) can  be-  ca  le- ul  a t e-il  elire-ctlv  from  this  i n feirma  t ion  , since-  p('  1 

\ V X y 


13 


is  tlie  probability  of  the  joint  occurrence  of  the  slope  pair  ,,  . 1 he 

11  X y 

Tormula  for  calculating  a joint  density  function  is; 


P(x^.yj) 

y>('  . ) = I T 

-X  y |j(x^,y^)l 


p(x  ,y  ) 
n 11 


J(x  ,v  ) 
n n 


where  the  Jacobian,  J,  is 


J(x,y)  - 


3 K ( X , y ) 
ox 

0 h ( X , y ) 
ox 


a ^ ( X , y ) 

ay 

o li  ( X , y ) 

oy 


CM) 


(35) 


and  tile  (x  ,v  ) ...  , (x  ,v  ) are  all  real  solutions  of  C =0  and 

1 ’ ■ 1 n ti  “X 

= O.  If  we  let  fj(x,y)  - 1 and  h(x,y)  = then  substituting  in 

the  .Jacobian  we  {lat  : 


J ( X , y ) 


‘ XX  ^xy 
^ xy  ' yy 


( .'it; ) 


•Since  the  perturbi'd  reflecting  surface  is  a surface  of  revolution, 
accordin';  to  .our  assumptions,  w<'  chanpie  the  variables  to  polai'  coordinates 
The  de  1' i va  t i ve.s  in  Cartesian  coordinati's  in  terms  ot  de  I’i  vm  t i ves  in  polar 


coordinates  are: 


Tile  .I;icobi:in  tuvns  out  to  l)c  iiulepenclont  of  c 


:i  nd  t Ini' 


j = U,.,,  '->■  'I'l 


p'  ..  > 

\ \ 


’’ll 


p { X , V ) 
n n 


(f  ) ^ (r  ) r 
11  r n n 


whore  r 


p(>:  ,v  ) = 1 nroii  = !/(•  i!  )‘ 
i ■ i max 


r a ff'  solutions  of  1 (r)  o. 

n “ r 


liocauso  of  the  nature  of  t )ie  per  t urlia  t ion  , we  find  t lui  t each  slope 

occurs  at  most  twice;  tlierelorc*,  n = 2.  We  also  notice  that  a maximum 

slope  always  occurs  lietween  tlie  two  solutions.  'the  radius  of  this 

maximum  slope  is  found  prior  to  any  calculation.  Tlie  solutions  of 

' (r)  are  then  localized,  and  r mav  fie  found  tiv  iteration.  'flic  function 

r 

' (r)  is  a linear  interpolation  on  data  points,  and  ' (r)  is  found 

>■  • rr 

analvtically  from  t fie  interpolation  formula. 


It  is  important  to  fiave  a ^ood  metfiod  of  interpolation  for  tlie  re- 
flectini;  surfac<‘  fiecause  we  are  not  (fea  1 i ny  directly  w 1 1 li  t fie  surface 
fiut  w i t li  its  slope,  1 fi  i s (|uaiitity  is  seiisitii'e  to  minor  cfianye^  in  t fie 
metfiod  of  i n t e rpo  1 a t 1 ny  fxtv.eeii  points.  'file  metfiod  adopted  i'  re  1 a 1 i vi' 1 y 
‘-impfe  and  introduces  \-erv  fittle  interpolation  error.  It  i a linear 

i n I e r [io  1 .1 1 i on  on  tlie  di  I’ivatiie  of  tlie  surface  w i t li  respect  to  t he 
logarithm  ol  tlie  di-lance.  tlie  derivative  w i t li  n-speet  to  loyaritlimic 
di-t.inc<'  i-  fouri'l  dirietl\  from  file  data,  wfiicli  consi.-ls  ol  tallies  of 
ri'lleetiny  -url.ici  lieiylit,  ",  a-  ,i  function  of  di-tanee  I roi-i  tlie  ci-ntii 
o I t 111  hur  - t , r . I li.i  \ I 
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I.ineor  interpolation  usint;  lo^nrithmic  distanci'  i^  ' riormod  on 
these  smoothed  data.  Thus 
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['or  Iom;  r^  loK  r ' Iok  r.  i = J,  n,  and 
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Ihe  set  of  sloiics  determined  tliis  way  is  inteprated  to  eliec-k  tlie 
apreement  with  t hc’  oripinal  rellectinp  surlaee.  It  t lie  smootliinp;  is 
too  severe,  additional  data  points  may  be  adcled  . 

A fCw  illustrative  calculations  ha\’e  been  made  usinp  l<|.  (dl),  with 

V chosen  to  pive  1 kU'  power  radiated  into  a liemisphere  ov<  r perlectlv 
o 

conductinp  pround,  lor  numbers  o!  d . a ~km  altitude  d-vl!  de  t ona  t i o.i  s , at 

-Cl  -li 

an  averape  diiisity  ol  10  km  , in  the  vicinity  ol  a 1 ooo-km  daytime 
path.  The  ell'ects  on  the  |iropapation  enviroiimetit  ol  a s 1 np  1 e such  di  t ona 
t ion  lU'fe  describid  in  detail  in  Sections  1 1 -A  and  111.  The  results  ol 
the  i neoherent -1  ie  1 il  calculations  are  summ;ir-i/ed  in  Tabli  b,  I rom  '^hich 
t.he  root -mean-squa  ri'  ( rms  ) maptiitude  ol  tlv  i iKuihi-rent  1 y scattered  liild 
, seen  to  be  in  the  dO  to  lo  V m ratipe.  Ihe  lields  luive  been  doabli'd 


from  the  values  tjiv'C'n  directly  by  Kq . CU)  to  account  for  tlic  assumed 
ground  cha rac tori st ics . 

■fal)le  (> 

Magnitude  of  ItMS  Vertical  Component  of  Incoherently 
Scattered  I-’ield  at  lOOO-km  Distance  from  l-k\V 
Source  for  a Detonation  Density  of  10“*’  km“- 
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I3y  comparison,  the  ambient  (colierent)  field  for  one-liop,  sky-\va\e 
propagation  over  sea  (a  highly  conducting  ground)  to  the  same  distance 
was  calculated*  to  l)e  about  OHO  . V/m  at  20  kll7 . Since  the  colu'rent  field 
is  essentially  undisturfied  at  this  detonation  density  (cf.  I'able  11),  tfiis 
va  1 U('--wh  icti  is  a pproxi  ma  1 <' 1 y 20  dl!  greatei-  tiutn  t lu’  lai-gest  (2-min)  va  1 lu 
for  the  rms  incohei'ent  field  at  Utis  frequency  — is  also  clia  rac  te  r i s t ic  of 
the  coherent  field  in  the  disturbc'd  environment.  .As  wo  noted  above, 
ttiis  smallness  of  tlie  incoherc'nt  field  I'elative  to  tli(>  coherent  field 
is  only  to  l)c>  expected  at  sucli  low  dt'tonation  (kuisities.  fhis  result 
should  not,  tiowever,  fie  taken  to  ctui  rac  tc  ri  x.<'  all  cases. 


Tile  ambient  field  stixuigth  was  determinerl  using  the  reflection  coeffi- 
cient given  by  l.(|  . (2S)  for  t lie  /ero-deiisitv  angle  coelficient  from 

lafile  .'i , in  coniunction  witfi  1 lie  wave-hop  pa  t li- i n t eg  ra  1 \-aIues  otitained 
tiy  Berry  and  ChrismaiT  for  a liO-km  r<f  feet  ion  height. 
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CONCLUSIONS  AND  KKCOMMKNDAT IONS 


'I'lio  calculations  completofi  here  indicate  that  hit;h-yicld  (MT  raiific), 
low-altitude  (lew  km)  nuclear  detonations  at  densities  of  10  km  “ or 
greater  have  substantial  effects  on  VI,1'/LF  propagation.  Such  densities 
of  detonations,  alttioufjti  hitch,  arc  not  unrealistic  for  lartcc-scale  attack 
scena  r ios . 

1 he  maior  effect  found  at  hifcli  d<?tonation  densities  is  an  increase 
in  t lie  rate  of  attenuation  of  tlic  coherent  ly -re  fleeted  propa  tea  t i nt;  wave, 
relative  to  that  in  the  undistui'bed  environment,  as  a result  of  the 
ionospheric  boundary  routchness  produced  by  debris  radiations.  The  consc- 
ciuent  attenuation  rates  are  fcrcatest  in  the  disturbed  nitcht  environment, 
as  a result  of  the  larjce  heifc'ht  ranjee  over  wliicli  reflection  occurs  for 
this  environment.  In  contrast  to  the  effects  produced  by  h iccfi-a  1 1 i t ude 
detonations,  daytime  attenuation  rates  in  the  disturlied  environment  a I’e 
s i ten  i f ican  t ly  less  than  those  at  ni<cht  . 

A method  has  been  presented  foi-  determination  of  the  efICcts  of  many 
randomly  disturbed  detonations  in  terms  of  t iiosc'  produced  by  a sintcle 
detonation.  This  method  assunu's  that  the  eflT'cls  at  all  but  a small  (and 
to  lirst  order  iu'k  1 i tcib Ic- ) traction  of  the  alfecte'd  retcion  are  determined 
primarily  by  the  nearest  ele  t eina  t i ein  tei  euich  peiint  ei  f the  reicion.  Ceimbi- 
natiein  e>f  this  assumptiein  with  a statisticeil  eiescri  pt  iein  of  the'  eletonatiem 
leieatieins  icre'iitly  simplifie'S  the  eissesswe'iil  eif  the'  I'livi  reinnu'Ht  . 'this 
appreiach  shenilel  alsei  be'  applicedilc'  t ei  eitlier  type  s of  nucle'ai'  e'ffe'e'ts  in 
aelelition  tei  that  ceins  1 ele' reel  here'. 

A feicteir  neit  e'\;im  i ne'el  epia  n t i I a t i ve' 1 y at  s i un  1 I i e-ant  tie' t ona  t i ein 
ele'iisi  tie's  is  the'  fi'eiclioii  eil  the-  e't  ihe  ■ I'e  ■ n t -s  i yna  1 leiss  t h;i  t peies  irtei 
the  i nceilu' ren  t 1 V se';i  t t e ■ re  el  wiive'.  ’fh  i s fr.ietiein  nuiy  neit  represe  ii  t a re-a  1 
lee^s  tei  the'  e'timmun  i e'a  t i ein  s syste'iii  unele'r  seime'  e' i re'umst  a nce'S  . More'  siteni- 
'irantlv.  the  i ne'ohe  ren  t 1 v se'a  t t efeel  wtive'  weiulel  be'  e-xpe-cte'el  •ceiie'rallv  tei 


I'cmain  effective  as  a .lamniinK  siKnal.  Clearly,  the  distribution  of  loss 
tietween  Ijoumlary  at)soi'ption  and  incotierent  scatter,  ami  the  utility  of 
the  latter  as  desired  and  as  .jamming  signal,  merits  further  study. 

The  calculations  reported  liei-e  must  be  considered  preliminai'y  in 
several  oilier  particulars.  Multiple  scatteriiif^  and  shadowing  have  not 
tieen  taken  into  account.  Further,  the  accuracy  of  the  statistical  model 
(wliicli  expix'sses  tlie  effects  of  tlie  ensemtile  of  detonations  in  terms  of 
tliose  of  a siiif^le  detonation)  has  not  lieen  determined.  These  aspects  of 
tlie  modeling-  sliould  all  lie  examined  and  suitaiile  laiio '’a t ion  of  tlie 
models  undertaken  if  necessary. 

I inally,  consideration  sliould  be  jiiven  to  incorporation  of  these 
models,  at  a suitable  level  of  detail  and  accuracy,  into  the  1)N.\  nuclear- 
effi'cts  codes  for  treating  VbF/bF  propaf^a  t ion . 
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